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Recent ARPES measurement observed a large afc-axis gap anisotropy, A(0, 7r)/A(7r, 0) = 1.5, in 
clean YBa2Cu307_<s. This indicates that some sub-dominant component may exist in the d x 2_ y 2- 
wave dominant gap. We propose that the interlayer pairing tunneling contribution can be determined 
through the investigation of the order parameter anisotropy. Their potentially observable features 
in transport and spin dynamics are also studied. 



PACS numbers: 74.20.-z, 74.72.Bk, 74.25.Ld 

The symmetry of the order parameter continues to be 
a fundamental issue in the studies of high-T c cuprate su- 
perconductors. While the order parameter of cuprates 
is extensively believed to be dominated by d x 2_ y 2-wave, 
whether there exists a sub-dominant component remains 
an open question [Q. Recent angle-resolved photoe- 
mission spectroscopy (ARPES) measurement H has re- 
vealed a significant superconducting gap anisotropy in 
clean untwinned YB^CuaOy-^ (YBCO). The excitation 
gap along the k y axis is found to be 50% more than 
that along the k x axis. This strongly suggests in YBCO 
that some sub-dominant component is involved in the 
d x 2 _j,2-wave dominant order parameter. It is important 
to ask what the symmetry and the mechanism of the 
sub-dominant component are. 

In terms of structure and electronic properties, 
YBCO differs from La 2 - K Sr 2 ,Cu0 4 + < 5 (LSCO) or 
Bi2Sr2CaCu208+,5 (BSCCO) in one crucial aspect. Due 
to the existence of one-dimensional (ID) CuO chain, the 
system is orthorhombic rather than tetragonal. The per- 
fect £>4/j symmetry is broken which accounts naturally 
for the additional symmetry mixed in the order param- 
eter. In this case, the most favorable order parameter 
is the so-called d+s model, Ak ~ cos(20) + s [0 is the 
azimuthal angle of k on the Fermi surface (FS)], which 
arises providing that the in-plane pairing interaction is 
given by a separable g(k, k') = — Vf(k)f(k') with V > 
and /(k) = cos(2</>) + s. The nodal lines of this gap are 
deviated from the usual diagonals (k x = ±k y ) and conse- 
quently the gap magnitude is anisotropic between the k x 
and k y axes. The ARPES data, A(0, ir)/A(n, 0) = 1.5, 
would correspond to s = —0.2 in this simple d+s model. 

Another complication occurs in YBCO (and also in 
BSCCO) because it involves more than one conduct- 
ing layer within a unit cell. This enables the interlayer 
Cooper pair tunneling in the superconducting state. Fol- 
lowing Chakravarty et al. || , a superconducting bilayer is 
considered to which an interlayer pair tunneling Hamil- 
tonian 

H± = ^2 Tj(k)a\ k1 -a\_ kl a 2 -kia2ki + H.c. (1) 

k 

coexists with the usual intralayer pairing interaction 
term. Here the index 1 (2) denotes for layer 1 (2) and 



Tj(k) is the interlayer Josephson pair tunneling cou- 
pling. The momentum-conserved interlayer pair tunnel- 
ing term (|l|) can be generated from the interlayer single- 
particle hopping ~ ^ k tj_(k)aJ k a2kTi where the hop- 
ping integral i_L(k) can be determined from the band 
structure M. Thus Tj(k) = t^(k) 2 /t is usually as- 
sumed with t the nearest-neighbor hopping. As pointed 
out by Chakravarty and Anderson J6|, due to the highly 
non- Fermi liquid nature of cuprates, the coherent inter- 
layer single-particle hopping is blocked, which neverthe- 
less gives rise to the interlayer pair tunneling. 

Combining ([!]) and the intralayer pairing term, one can 
write down the BCS gap equation for layer i (1 or 2) 

A k i = -^.g(k,k')(ai_ ka a ik , T ) + 7j(k)(a.,_ ki a., kT ), (2) 

k' 

where j ^ i and the angular bracket denotes the anoma- 
lous mean-field average. It should be stressed that 
the two (Cu02 plane) layers considered here for YBCO 
are mediated and renormalized by the third CuO chain 
layer and appear to be orthorhombic instead of tetrago- 
nal. The issue is then how the interlayer pair tunneling 
plays the role between the two orthorhombic Cu02 lay- 
ers. For two identical layers, we assume (ai-k|aikt) = 
(a 2 _kia2kT) = (A k /2.E k )tanh(£: k /2T) by symmetry, 
where £" k = [£ k + A k ] x / 2 is the quasiparticle excitation 
spectrum with £ k the particle band energy and A k the 
overall gap self-consistently determined by (||). At T = 0, 
Eq. (U) is reduced to (layer index is redundant) 

A k = A /(k) + A k ^, (3) 

where Ao = VY^k" A k /(k)/2_E k with lu c an appropriate 
BCS cutoff frequency. At or near the FS (|k| ~ kp), 
Ek ~ |A k | and (||) is further reduced to 

A k = A /(k) + ^sgn[A k ]. (4) 

Due to the term sgn[A k ] , it sets a constraint in (Q) that 
physical solutions of A k arise only when Tj (k) preserves 
the symmetry of /(k) which in turn is determined by 
the nature of the layer pairing interaction. When the 
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FIG. 1. The gap magnitude as a function of the azimuthal 
angle on the FS for a superconducting orthorhombic bilayer 
with interlayer pair tunneling [using (El)]. 



symmetry of Tj(k) is different from that of /(k), the 
overall gap could exhibit a discontinuous phase on the FS. 
This would imply a physically unstable state. To fulfill 
the symmetry requirement, we thus argue that ij_(k) oc 
[cos(2</>) + s] 2 andTj(k) oc [cos(2 cj>) + s] 4 assuming /(k) = 
cos(20) + s for the present superconducting orthorhombic 
bilayer. More precisely, the overall gap function in (Q) is 
given by 



A k = A [cos(20) 



-y [cos(2</>) + s] 4 sgn[cos(2</>) + s] 



(5) 



where Tj denotes the strength of the interlayer pair tun- 
neling. Eq. (|) shows that not only the orthorhombicity 
s, but also the pair tunneling strength Tj can be deter- 
mined through the investigation of the gap function. It is 
worth emphasizing that when s is small, the a6-axis gap 
anisotropy could still be large, as long as Tj is significant. 

The band structure calculation and experimental 
evidences j^H usually assume t±(k) cx (cos k x — cosk y ) 2 
[or oc cos 2 (2(f)) in the continuum limit] for high-T c 
cuprates in a square lattice of perfectly tetragonal sym- 
metry. This justifies our proposal ^i(k) oc [cos(2<^>) + s] 2 
in the tetragonal limit (s = 0). 

Fig. [l] shows the gap magnitude as a function of 4> 
for two extreme cases: s — — 0.2, Tj = and s = 
-0.06, Tj = 7.6A . (In fact, A is a function of Tj 
through the self-consistent gap equation.) These two 
cases all give A(90°)/A(0°) = 1.5 in regard to the 
ARPES result. Several important features are noted for 
the anisotropic gap. The nodal angle, deviated from the 
diagonals, is now determined by 



6 nodc = i arccos(-s) ~ ^ 



(for small s) 



(6) 



in the first quadrant of the FS. Measurement of the nodal 
angle thus reveals the value of s. With the knowledge of 



s, one can then compare the gap magnitude between the 
two antinodes 



|A(90°)| (1 - s)[l + r(l - sf 



|A(0°)| (l + s)[l + r(l + s) s 



Tj 
2A : 



(7) 



to obtain the value of Tj (in unit of A ). Alternatively 
one can also study the gap slope near the nodes 



d[A(0)/A(O°)] 



= 2 



(l + s)[l + r(l + s) 3 ] 



(8) 



which is 2 for a pure d 2 ,2_ ?/ 2-wave gap (s = Tj = 0). 

If the instrumental resolution is fine enough, ARPES 
would be the best probe to explore Tj and s through 
the direct measurement of gap anisotropy. In the fol- 
lowing, we study longitudinal ultrasonic attenuation and 
inelastic neutron scattering spectra which are also useful 
probes on these issues. 

Longitudinal Ultrasonic Attenuation — Ultrasonic at- 
tenuation is a directional probe and thus powerful to 
study the gap anisotropy. It has been successfully used 
to study the order parameter symmetry in superconduct- 
ing heavy-fermion UPt 3 @ and high-T c analog Sr 2 Ru0 4 
compounds 113]. In clean YBCO where the scattering 
is in the ballistic limit H], the longitudinal ultrasound 
attenuation in the superconducting state is proportional 

to ra,Ei 



a s (q,T) oc 



k 



df(E k ) 



i 6 



ak 



A, 



dE k 
: dk 



f 2 



(9) 



where q is the wavevector of the propagating phonon and 
/ is the Fermi distribution function. We shall calculate 
@ for one single-layer with the input of the interlayer 
pair tunneling gap in (|5|). This is sufficient if no other 
coupling or vertex correction is considered within the two 
layers. Inspection of (j^ shows that a s is governed by a 
delta function, weighted by the FS sum according to how 
small the gap is. For an isotropic s-wave superconductor, 
the second term in the delta function vanishes, i.e., only 
the portion of the FS perpendicular to q contributes. For 
unconventional superconductors such as YBCO, the sec- 
ond term does contribute but is typically small (oc A/cf 
with A the maximum gap and cf the Fermi energy). By 
ignoring the second term in the delta function, Eq. (JsJ) is 
simplified to 



a s _ (2J(k F -q)/(|A k |)) 
a n (<5(k F -q)) 



(10) 



between the superconducting and normal states, where 
the angular bracket denotes the FS average. Since only 
two points on the (two-dimensional) FS satisfy k^-q = 0, 
Eq. © is simply solved to be a s /a n = 2/(|A(0)|) = 
2/[exp(|A(0)|/T) + 1], where is related to (f> q (the angle 
of q) by - <j> q = it/2. 
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FIG. 2. The ratio of the longitudinal ultrasonic attenua- 
tion in superconducting to normal state as a function of the 
angle of sound wavevector. Top (bottom) frame corresponds 
to s = -0.2 and Tj = (s = -0.06 and Tj = 7.6A ). The 
curves from top to bottom are for T/T c — 0.5, 0.3 and 0.1. 



In Fig. a s /a n vs. <j> q are shown at different tem- 
peratures for those two cases in Fig. We have as- 
sumed 2A(0°) /T c — 3.52 and the gap magnitude remains 
roughly the same for T < 0.5T C . [Of course, the exact T 
dependence of gap can be worked out through (j^), which 
is not done in this paper.] The maximum a s /a n appear 



at a specific angle <f> q at 



7T 

2 



(11) 



where the nodal angle </> n odc is given in (g) and we limit 
< </>max, ^node < 7r/2. Since </> n odc depends on the value 
of s, determination of </> max in ct s /a n thus determines s. 
In addition, the ratio at <f> q — to 7r/2, 

a s /a n {(j) q = 0) 



a s /a n (4> q = 7r/2) 

1 + exp [jf ((1 + s)[l + r(l + s) 3 ])] 
l + exp[^((l-s)[l + r(l-s) 3 ])j 



(12) 



allows one to extract the strength of Tj. Furthermore, 
near the critical value, the slope 



d(a s /a n ) 



d<t> q 



A 



s/T 



T (l + s)[l + r(l + s) E 



(13) 



which gives one more way to look into Tj and s, as well 
as Aq. 
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FIG. 3. Inelastic spin excitation spectra with the momen- 
tum scan along the direction shown in the inset. Top (bot- 
tom) frame corresponds to s = —0.2 and Tj — (s = —0.06 
and Tj — 7.6Ao). The curves from top to bottom are for 
cj/A(0°) = 1.0,0.8 and 0.6. 



Spin Excitation Spectra — Recent inelastic neutron 
scattering (INS) experiments have reported many inter- 
esting results for high-T c cuprates regarding the commen- 
surate and incommensurate (IC) peaks at or near the an- 
tifcrromagnetic (AF) wavevector Qaf = (7T, 7r ) 0. For 
YBCO or BSCCO, the appearance of the commensurate 
peaks are likely due to some kind of resonance, while the 
IC peaks are associated with the dynamical local nesting 
effect of the band structures ||T5|— [l7|] . 

Theoretically the INS spectra is proportional to the 
imaginary part of the spin susceptibility and, for sim- 
plicity, we study the irreducible BCS spin susceptibility 
Xg(q, cj) for one single-layer. In Fig. [| Imxo are calcu- 
lated at T — for some fixed frequency to and momen- 
tum q scanned along the direction shown in the inset. 
The primary interest is to see how the gap (|5|) affects 
the INS spectra. To study the IC peaks, a tight-binding 
band £k = — 2t(cosk x + cosk y ) — At' cos k x cos k y — fi 
is inevitably used, where t and t' are respectively the 
nearest-neighbor and next-nearest-neighbor hopping and 
fi is the chemical potential. Typical values of t' — — 0.25i 
and [i = — 0.65i are employed. For the gap, cos(20) in 
(||) is replaced by (cosfcj, — cos/cj;)/2 in the present lat- 
tice case and we assume the gap magnitude at fc^-axis, 
A(0°) = 0.3i. The same two cases in Figs. |l| and || are 
studied. For the Tj = (larger s) case, the nodes are 
shifted away from the diagonals and consequently the 
nesting effect is highly anisotropic: the IC peak parallel 



3 



to q y axis (along the q scan route chosen) is more intense 
than the one parallel to q x axis. At lower frequency, 
the latter could completely disappear before the nesting 
effect can set in. For the strong Tj (small s) case, in 
contrast, the nodes remain closer to the diagonals and, 
as a result, the nesting effect is roughly isotropic. The 
1C peaks are nearly symmetric along the q route, regard- 
less the change of frequency. These are the key features 
which can be used to distinguish the strong and weak Tj 
(or small and large s) cases. 

When the commensurate peaks are also considered, the 
Random phase approximation (RPA) corrected spin sus- 
ceptibility x s (q.w) = xo(q>^)/[ 1 -^(q)x§(q,w)] is often 

studied, where V(q) is usually modeled by an "AF" in- 
teraction V(q) = —J(cosq x + cosg y )/2 (J > 0). Since 
here we are only interested in the IC peaks which appear 
beyond the AF resonant regime, it's adequate to study 
Xq. The line shape in Imx s (q, lu) will be very similar to 
those in Xp(q, lu), apart from different intensity. 

Mook et al. Jlq] have recently reported the ID nature 
for the IC peaks in underdoped YBa2Cu30g.6- The IC 
peak intensity at (ir — S, n) is found to be stronger than 
that at (tt, it — 6). While the results are referred to favor 
the formation of the dynamical stripes for that particular 
doping [p| , an alternative explanation to these data is to 
take into account the in-plane gap anisotropy associated 
with a homogeneous Fermi liquid Jh| . The data of Mook 
et al. seems indicating an intermediate value of Tj (and 
an intermediate and positive s) on YBa2Cu30g.6, when 
comparing to the results in Fig. ||. It is important to 
have INS performed at different frequencies in order to 
clarify the gap anisotropy. 

In summary, we propose that the long-thought inter- 
layer pair tunneling effect can be explored through the 
detailed studies of the gap anisotropy in YBa2Cu307_a. 
Considering an orthorhombic superconducting bilayer 
with interlayer pair tunneling, one can simultaneously 
determine the in-plane orthorhombicity and the strength 
of the pair tunneling. For probable probes, we study 
longitudinal ultrasonic attenuation and inelastic neutron 
scattering. 

Finally, we comment on the effect of pair tunneling 
on the collective modes. For a superconducting bilayer 
coupled by the pair tunneling, apart from the gap renor- 
malization, the system will exhibit the characteristic 
in-phase and out-of-phase phase modes associated with 
the order parameters on the two layers. The in-phase 
Anderson-Bogoliubov phase mode has a phonon disper- 
sion, lu 2 = v 2 F q 2 /2, which is lifted to plasmon mode 
when long-range Coulomb potential is included. These 
in-phase modes, which are intrinsic characteristics of a 
superconductor, are independent of Tj. In contrast, the 
out-of-phase phase mode has an optical phonon disper- 



sion LU 



■ w 



Vpq /2, where luq depends on the relative 
size of planar pair interaction V and Tj [ p(i[ . It is also 
interesting to study Tj through the search of the out-of- 
phase phase modes, such as those recently observed in 
SmLao.8Sro.2Cu04_5 with two different Josephson cou- 



plings ||. 
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